STAR 



ADAPTION OP A CORRECTOR MODULE 
TO THE IMP DYNAMICS PROGRAM 


Avco Systems Division 

201 Lowell Street 

Wilmington, Massachusetts 01887 


N73-10238^\ 


Onclas 

46607 


( NAS A-CR-1 30108) 
MODULE TO THE IMP 
Report and user's 
Wilmington, Mass.) 


ADAPTION OF A CORRECTOR 
DYNAMICS PROGRAM Final 
manual (Avco Corp. , 

Sep. 1972 103 p C5CL 

09B 


13/08 


September 1972 

Final Report and User's Manual 
Contract NAS5- 11881 
AVSD-0335-72-RR 


Prepared for 

GODDARD SPACE FLIGHT CENTER 
GREENEELT, MARYLAND 20771 




ADAPTION OF A CORRECTOR MODULE 
TO THE IMP DYNAMICS PROGRAM 


Avco Systems Division 

201 Lowell Street 

Wilmington, Massachusetts 01887 


September 1972 

Final Report and User's Manual 
Contract NAS5-11881 
AVSD-0335-72-RR 


Prepared for 

GODDARD SPACE FLIGHT CENTER 
GREENBELT, MARYLAND 20771 


* 

I 



TABLE OF CONTENTS 


Page 


1.0 INTRODUCTION 


1 


2.0 TECHNICAL DISCUSSION 5 

2.1 Performance Criterion 8 

2.2 Search Logic 10 

3.0 PROGRAM DOCUMENTATION 16 


3.1 Subroutine Descriptions 17 

3.2 Definition of Labeled Common 34 

3.3 Definition of Input 65 

3.4 Description of Output 80 

3.5 Overlay Structure 8l 

3.6 User's Instructions 86 


4.0 DETERMINATION OF ANTENNA POSITIONS FROM ATTITUDE DATA 90 


4.1 Simplified Equations of Motion 91 

4.2 Illustrative Corrector Runs 93 

4.3 Conclusions 96 


5.0 REFERENCES 


98 


-ii- 



LIST OF FIGURES 


Figure No . Page 

1 Corrector Schematic 3 

2 Typical Search Sequence for Two Parameter 

Case 13 

3 Search Completion Criterion 14 

4 MAIN Block Diagram with Corrector Added .... 18 

5 Labeled Common - SUBROUTINE Interaction .... 35 

6 Overlay Tree Structure 83 

7 Geometry of X Configuration 91 

8 Deformations for Generalized Coordinates ... 93 


-iii- 



1 . 0 INTRODUCTION 
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The corrector module of HAEIOS Program (Reference 1) and 
the IMP Dynamics Computer Program (Reference 2) have been com- 
bined to achieve a data fitting capability with the more general 
spacecraft dynamics models of the IMP Program. The IMP Dynamics 
Program models spacecraft dynamics for satellites having long 
flexible booms . The resulting equations of motion are a system 
of ordinary and partial differential equations. The use of mode 
shapes to represent the deformations of the flexible booms con- 
verts the partial differential equations into ordinary differen- 
tial equations. The resulting global system of ordinary differen- 
tial equations is nonlinear, time-dependent in some of its forcing 
functions and, at least piecewise, time-invariant in its system 
parameters such as antenna boom stiffness, damping, thermal prop- 
erties, etc. 

In order to start the IMP Dynamics Simulator > initial values 
of state variables are required. These state variables include 
central core attitude angles and rates, antenna tip displacements, 
and velocities, etc. Estimates are also required of the system 
parameters. The closeness of fit of the Dynamics Simulator depends, 
among other things, on the accuracy of these estimates. 

The general corrector scheme devised to supply "best" esti- 
mates of the uncertain system parameters and missing initial con- 
ditions is an iterative procedure. The overall logic, shown in 
Figure 1, is applicable to practically all corrector or optimiza- 
tion schemes. The comparator computes the difference between 
simulated behavior and measured data in order to generate error 
signals . The performance criterion combines all error signals 
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into a single scalar variable which expresses the quality of the 
fit between simulated and measured data. The Corrector allows a 
choice of different structures for the performance criterion, 
i.e., weighted or unweighted and absolute values of errors or 
squares of errors. The weighted, squared-error performance 
criterion structure can be thought as (and is in fact equivalent 
to) the introduction of a diagonalized error covariance matrix, 
which is often used in filtering problems . 


A Priori Estimates 



Figure 1 - Corrector Schematic 
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The remainder of this report describes the properties of 
the corrector. The technical discussion is concerned with the 


description of the performance criteria and the search logic for 
parameter estimation. The final section gives a description of 
the modifications made to add the corrector to the IMP Dynamics 
Program. This section Includes subroutine descriptions, common 
definitions, definition of input, and a description of output. 
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2.0 TECHNICAL DISCUSSION 
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A variety of search algorithms exist and were examined before 
selection of the method used in the Corrector. The Rosenbrock 
rotating coordinate scheme which belongs to the so-called direct 
search class was selected for the following reasons: 

1. Because of the complex nonlinear and time dependent 
nature of the Dynamics Simulator, forming mathemati- 
cal expressions for the derivatives of the state 
equations needed for gradient search methods is a 
very cumbersome task. Of course, it is possible to 
estimate numerically such derivatives, but such 
estimates are very sensitive to noise amplification 
and, moreover, would require a great number of simu- 
lation runs for a high dimensional state-space. 

Direct methods avoid exact estimates of gradients. 

2. Direct search methods, especially of the rotating 
coordinates variety, are quite independent of the 
nature of the criterion surface in parameter space, 
i.e., they are applicable to a great variety of 
criterion surfaces that do not even have to be con- 
tinuous. While gradient techniques may be faster 
in convergence on nearly linear or nearly quadratic 
portions of surfaces (as may be the case near the 
optimum), they may be much slower on more irregular 
surfaces or even diverge completely. The Rosenbrock 
method, on the other hand, has forced improvements 
at nearly every step. Because of the complexity of 
many satellites, it is not, a priori, evident what 
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are the characteristics of their criterion surfaces. 
Therefore, the Rosenbrock method, with its wider 
field of applicability, appeared to>'be a safer 
choice . 

.3. The direct search algorithm adapts itself very 

easily to changes in the Dynamics Simulator and to 
changes in search parameters, even if these changes 
are not trivial. 

The Rosenbrock search technique, as described in Reference 3 , 
searches for the greatest or least value of a function of several 
variables. The method is an algorithm for selecting trial values 
for the input parameters of a given system model in such a way 
that some error function of the performance of the system tends 
to be optimized. The technique has been referred to as the 
method of "rotating coordinates" because of the manner in which 
the direction of the variables are perturbed after a series of 
successful and unsuccessful searches have been made. Rather than 
varying the input parameters of the system model one at a time, 
the method rotates the search coordinate system, which initially 
coincides with the parameter axes in parameter space, so that 
one search axis points in the most favorable direction as deter- 
mined from the previous search. The remaining search axes, which 
are mutually orthogonal, are obtained with a Gram-Schmidt ortho- 
gonalization procedure. The search is continued along each of 
these directions one at a time using the logic which is described 
below, and then a new set of axes are developed. 


- 7 - 



The following list of initial conditions and spacecraft 
properties can be used as parameters in the optimization pro- 

. V.'/, 

cedure: 

1. Initial conditions on attitude and attitude rates 

2. Initial conditions on antenna displacements and velocities 

3. Initial conditions on damper angle and angular rate 

4. Antenna frame Euler angles 

3. Antenna modulus of elasticity 

6. Antenna temperature gradients 

7. Antenna lengths 

8. Antenna no load offsets 

9. Antenna damping coefficients 

10. Libration damper magnetic hysteresis decay factor 

11. Libration damper magnetic hysteresis saturation torque 

12. Libration damper magnetic hysteresis initial magnetic 
moment 

13 . Libration damper suspension spring constants 

14. Libration damper suspension stop angle 

15- Libration damper suspension viscous damping coefficient. 

2 .1 Performance Criterion 

An error function or performance criterion is determined 
from a comparison of actual flight data and simulated data over 
a specified time span. The flight data is read from a tape and 
stored in common for use in evaluating the performance criterion. 
Plight data must be supplied at a constant time interval over 
the specified time span (for instance, one minute intervals). 

The flight data is stored in a vector having time and the data 
for each data channel at every equally spaced time point in the 
specified data span. 
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The performance criterion is constructed for each data channel 
during each trial simulation. At each data time a simulated data 
vector is formed. An error vector is computed as the difference 
between the flight data and the simulated data. Provision is 
made to weight and normalise the error vector to achieve a weighted 
and normalised error vector if desired, as follows. 



( 1 ) 


Where: 


Weighting factor for the Oth data channel. 

£ .= Difference between flight data and simulated 
6 data for the Och data channel at the h time. 

CL- •= Normalises to the corresponding flight data. 

V If the flight'data is less than the normalisation 
constant, Ac* SiJ- is normalised to A^. 

This avoids division by zero. 



"= / “y/ > A; 1 

Ac if I a,y\ Ai 


( 2 ) 


The accumulated error in each data channel is computed by adding 
either the absolute value or the square of to the appropriate 
data channel residual. 


i 

* ! 


Absolute Value Criterion 


Squared Criterion 


( 3 ) 
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At the end of the specified data span, the total performance 
criterion is obtained by summing the accumulated residual errors 
for each data channel. 


J= l R, 

L 

2.2 Search Logic 


(*) 


Each time the performance criterion is computed constitutes 
a "trial" . A trial is a "Success" if the performance criterion is 
equal to or less than the value obtained on the previous trial; 
otherwise it is a "failure" . Trials are made along one search 
axis in parameter space until there is a failure which has been 
preceded by at least one success . The length of step for these 
successive trials are determined as follows: 

i) An initial step size of e n is tried. The initial 

values of e n in each direction are input to the program; 
subsequent values are dependent on previous trials . 

In particular, if d n is the algebraic sum of all 
successful trials in the nth direction, then the 
initial trial of the next search in the nth direction 
will be ^ d 

ii) After a successful trial, the length of the previous 
step is multiplied by a constant oC (<* >i) and this is 
added to the previous value used. 


iii) If the initial step in a given direction is a failure, 
then the initial step is replaced by — /Sc,,^ (o< /3-c 
If this fails, the step is multiplied by an additional 
— etc. This procedure must result in«an eventual 
success since in the limit the search will return to 
the initial starting point. 


where o< < / 


After a set of trials has been completed on one search axis, 
the program searches along the next orthogonal search axis until 
all search axis have been treated. A new set of search axis is 
then calculated. All of the trials along the search axis and the 
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subsequent calculation of a new set of search axes is called a 
"stage" . The number of search axes correspond to the number of 
parameters being optimized. The rotating search axes are related 
to the parameter axes by the direction cosine matrix [o,J with 
which steps on a given search axis can be resolved into parameter 
changes. For the first stage [Cm ] is a unit matrix so that 
each step, e n corresponds to a change in only one of the system 
parameters. For each subsequent stage, a new direction cosine 
matrix is computed using the Gram-Schmidt procedure as follows: 


A 

c ° c ° 

Let 5 5 X • • • ' ' be the set of orthogonal unit 

vectors defining the directions in the original stage. Suppose 

that d n is the algebraic sum of all successful steps e n in 
i i 

the direction 5 , etc. Then define the set of vectors: 


A 

• s; 


K = 


etc. Then define 

the set 

of 

A 



+ • 1 • 

• • • <£» '■ 




A c 

<^4 % +• 1 ' - • ' 






The orthogonal unit vectors 



stage are now obtained using the following vector equations: 
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The new direction cosine matrix is obtained by taking the 

^ y 

inverse of the matrix comprised of the components of along 

the parameter axes . 

A typical search sequence for a two parameter case is shown 
in Figure 2 to illustrate the above procedure. Contours correspond 
to constant values of the performance criterion J. The first two 
stages shown correspond to 13 trials . 

In order to terminate the search, two controllable completion 
criteria are included. The first completion criterion terminates 
the optimization procedure if the performance criterion from 
three successive stages falls within a specified band. The second 
criterion terminates the procedure if the ratio of three successive 
performance criteria is less than a specified value, The basic 
features of the completion criterion are shown in Figure 3. The 
parameter D which defines the band for the first criterion may 
vary between zero and one. Specifying a value near 1.0 would 
emphasize the second criterion which is controlled by the value 
of the specified ratio, R. A value of D near 0.0 would make 
the termination dependent upon the magnitudes of successive J's 
as shown in Figure 3 . 

Specifying a small value of R near zero will terminate the 
optimization near a minimum of the performance criteria providing 
the convergence is not excessively slow. If the convergence is 
very slow, the first criterion will terminate the problem. 

The computer running time for the corrector is greatly 
dependent upon the input constants used for D and R of the search 
termination. The particular values used depends on the rate that 
the problem converges. Initially, the value of D should be small 
until the convergence rate can be observed. If the convergence 
rate is rapid, D can be set nearly equal to 1.0. If R ^ 0, 
the problem will terminate only after finding a minimum if D = 1.0. 
For initial studies, a value of R between 0 and 0.5 will terminate 
the search before a minimum was reached and conserve computer 
running time. The value of R should not be less than or equal to 
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UNIT VECTORS 
!•» STAGE 


UNIT VECTORS 
2 nd STAGE 



Typical Search Sequence for Two Parameter Case 




J(Y) 


J (Y + l) 


J(Y +2) 




DJ(Y) 

4 

< 

/ /n 

E-JIV) 

C 

f-r 

J(Y + I) 

1 

(l-D)J(Y) 




i 

i 


STAGE 

Y 


STAGE 

Y+l 


I 

STAGE 
Y + 2 


SEARCH IS COMPLETED IF BOTH OF THE FOLLOWING 
CONDITIONS ARE SATISFIED: 


CONDITION 1: (1 - D) J (Y) < J(Y + 1) 

(I - D) J (Y) < J (Y + 2) 

CONDITION 2: J(Y + 1) - J (Y + 2) < R 

J(Y) - J(Y + 1) 


Figure 3 - Search Completion Criterion 


-J( Y + l) 
J(Y+2) 
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zero since successive values of the performance criterion, J, 
are always less than preceeding values. 

As a further safeguard, the maximum number of trials conducted 
in a search is specified through input. Although the search may 
not be complete upon the execution of the specified trials, an 
indication of search progress can be obtained from the resultant 
printout. The user may then alter input and resubmit to obtain 
optimized parameters. The altered input will reflect the trends 
and results established by the incomplete search. 

The Rosenbrock method for finding the minimum or maximum 
value of a function suffers from the same limitations of most 
optimization procedures in that it assumes the function J (y) 
is unimodal, i.e., has only one minimum over the range of values 
of interest. In general, the program will converge to a minimum 
of J. One approach to exploring the uniqueness of this, minimum 
is to initiate the procedure with new initial values of the 
parameters; the process should always converge to the same optimum 
parameters if the function J (y) is unimodal. 
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3.0 PROGRAM DOCUMENTATION 
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3.1 Subroutine Descriptions 

The new subroutines are documented in the following 
paragraphs. The modifications made to existing routines are 
described below. The modified flow chart of the main routine is 
shown in Figure 4 . 

MAIN - The main routine was modified to control corrector 
operation . 

Cj^NVRT - This subroutine was modified to add a single 

coefficient for bending stiffness affecting all antennas. 

G0UT - This subroutine was modified to suppress printing 
of all corrector trials if desired. 

GPRINT - This subroutine was modified to suppress printing 
of all corrector trials if desired and to write 
test flight data sets . 

SJ^UT - This subroutine was modified to suppress printing 

of all corrector trins if desired. 

S0LAR - Modified to save space 

C0MBNZ - Modified to save space 
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ENTRY MAIN 








Figure 4 ~ MAIN Block Diagram With Corrector Added 
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SUBROUTINE NAME : ARFIC 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To preset control data when using standard 

flight data sets 


CALLING SEQUENCE: Call ARFIC 

LABELED COMMONS USED: ICNC0R, ICRPIN, IMAIN1, IRHIST, ISDARR, 

ISIMDT, IVARYP, XDBUG 


SUBROUTINES CALLED: No subroutines 


GLOSSARY OP VARIABLES: All variables defined in labeled common 



SUBROUTINE NAME: BL/CK4 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91.95 (OS) 

PURPOSE: To preset data 


CALLING SEQUENCE: Call BL/CK4 

r* 

LABELED COMMONS USED: PRCOM, IVARYP, EVARYP, BCNVRT, XDBUG , ISIMDT, 

IRHIST, IR0SIE, C0UFCN, CRjZfelE, ICNC0R, 

CROSY, TWRIT 


SUBROUTINES CALLED: No subroutines 


GLOSSARY OF VARIABLES: All variables defined in labeled common 
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SUBROUTINE NAME: CRINIT 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91.95 (OS) 

PURPOSE: This subroutine sets up the vector of parameters P(l) 

to be varied and the vector of initial variations E(l) 

CALLING SEQUENCE: Call CRINIT 

LABELED COMMONS USED: IVARYP, EVARYP, XINl, RP00L1, XIN2, CSTVAL, 

ANTENA, XIN3, CCNVRT, RTDIST, CSOLAR, CFNALP, 
CDMPER, RDAMPR, XDBUG, CFNDGB, ICRPIN, 

PARRA Y, BCNVRT, ICNC/R, IMAIN1, IP//L1 

SUBROUTINES CALLED: 

GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

X01(3 ) 



R*8 

Substitution symbol for 
PSI1, THET1, PHI1 

X02 ( 3 ) 



r*8 

Substitution symbol for 

/meg ( 3 ) 

X03(6) 



r*8 

Substitution symbol for 
ALFAE, BETAE, GAMAE , 

/mbc .( 3 ) 

X04(4) 



r*8 

Substitution symbol for 
BETLD, GAMLD, PHILD, 
DPHILD 

X05( 120) 



R*8 

Substitution symbol for 
A(10,3), AD/T(10,3), 
B(10,3), BD/t(10,3) 

X06 ( 120) 



r*8 

Substitution symbol for 
DIN(l-0,3), DIND/T(10,3), 
D0UT (10,3), D/UTDT (10,3) 

X07(30) 



r*8 

Substitution symbol for 
ALFAEK(IO'), BETAEK(IO), 
GAMAEK(IO) 


- 21 - 


SUBROUTINE NAME: CRINIT (Cone I'd) 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: 

CALLING SEQUENCE: 

LABELED COMMONS USED: 


SUBROUTINES CALLED: 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

X08(10) 



r*8 

Substitution symbol for 
EM0DLS ( 10 ) 

X09(l0) 



r*8 

Substitution symbol for 
TDIS(IO) 

X10(10) 



r*8 

Substitution symbol for 
ZLO(IO) 

xii(6o) 



r*8 

Substitution symbol for 
SK0A(1O,3), SK0B( 10, 3 ) 

X12(6) 



r*8 

Substitution symbol for 
DECAY, ZMD/, ZMDBJ0, ZK1D, 
ZK2D, PHIS 

X13 



r*8 

Substitution symbol for 
CNV 

Xl4 



r*8 

Substitution symbol for 
DT00 

X15 



r*8 

Substitution symbol for 
CDAMP(3, TO) 

Xl6 



r*8 

Substitution symbol for 
BSTIF 
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SUBROUTINE NAME: GRAM (C,D,NMAX) 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To calculate the transformation (C) which relates the 

'new' search axes to the original parameter axes 

CALLING SEQUENCE: Call GRAM (C,D,NMAX) 

LABELED COMMONS USED: No labeled common 


SUBROUTINES CALLED: No subroutines called 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

C 

1 1 

o 

M 

1 1 


R*8 

Transformation matrix to 
transform parameter 
changes on 'new' search 
axes to the original 
parameter axes 

D 

KJ 


R*8. 

Accumulated parameter 
changes on last completed 
stage 

NMAX 



1*4 

The number of parameters 
in the search vector 
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SUBROUTINE NAME: LDARAY 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To load the new parameter value into appropriate common 

array at the completion of a trial and 


CALLING SEQUENCE: Call LDARAY 

LABELED COMMONS USED: XIN1, RP//L1, XIN2, XDBUG, CSTVAL, ANTENA, 

XIN3, CCNVRT, RTDIST, CS/LAR, CFNALP, 
BCNVRT, CDMPER, RDAMPR, CFNDGB, ICRPIN, 

P ARRAY, ICNC/R 

SUBROUTINES CALLED: No subroutines 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

X01(3) 



R*8 

Substitution symbol for 
PSI1, THET1, PHI1 

X02 ( 3 ) 



r*8 

Substitution symbol for 
/MEG(3) 

X03(6) 



R*8 

Substitution symbol for 
ALFAE, BETAE, GAMAE, 

0mbc ( 3 ) 

X04(4) 



r*8 

Substitution symbol for 
BETLD, GAMLD, PHILD, 
DPHILD 

X05(120) 



r*8 

Substitution symbol for 
A ( 10, 3 ) , AD/T(10,3), 
B(10,3), BD/T(10,3) 

X06(120) 



R*8 

Substitution symbol for 
DIN( 10, 3 ) , DIND/T(10,3), 
D0UT( 10,3), 0/11101(10,3) 

X07(30) 



r*8 

Substitution symbol for 
ALFAEK(IO), BETAEK(IO), 
GAMAEK(IO) 
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SUBROUTINE NAME: LDARAY (Cone I'd) 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: 

CALLING SEQUENCE: 

LABELED COMMONS USED: 


SUBROUTINES CALLED: 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

X08( 10) 



R *8 

Substitution symbol 

em/dls(io) 

for 

X09(10) 



R *8 

Substitution symbol 
TDIS(IO) 

for 

X10(10) 



r *8 

Substitution symbol 
ZLO(IO) 

for 

xii(6o) 



r *8 

Substitution symbol 

sk/a(io,3), sk/b( 10 j 

for 

3) 

X12(6) 



r *8 

Substitution symbol 
DECAY, ZMD/, ZMDB/, 
ZK2D, PHIS 

for 

ZK1D, 

X13 



r *8 

Substitution symbol 
CNV 

for 

Xl4 



r *8 

Substitution symbol 
DT00 

for 

X15(30) 



r *8 

Substitution symbol 
CDAMP(3, 10) 

for 

Xl6 



r *8 

Substitution symbol 
BSTIF 

for 
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SUBROUTINE NAME: LSARAY 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To extract "simulator output from common and assign it 

to appropriate data channel 

CALLING SEQUENCE: Call LSARAY 

LABELED COMMONS USED: CANTNA, C0NSTS, ISDARR, RP00L1, SIMDAT, 

XDBUG, XIN1, XIN2 


SUBROUTINES CALLED: No subroutines called 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

X01(3) 



R*8 

Substitution symbol for 
PSI1, THET1, PHI1 

X02 ( 3 ) 



r*8 

Substitution symbol for 
0MEG(3 ), in degrees/sec 

XR2 ( 3 ) 



R*8. 

Substitution symbol for 
0MEG in rad/sec 

X03(6) 



r*8 

Substitution symbol for 
ALFAE, BETAE, GAMAE, 
0MBC ( 3 ) 

X04(2) 



r*8 

Substitution symbol for 
PHI, PHID in deg, deg/sec 

XR4(2) 



r*8 

Substitution symbol for 
PHI, PHID in rad, rad/sec 

X05( 120) 



r*8 

Substitution symbol for 
A(10,3), AD0T ( 10, 3 ) , 

B( 10, 3 ) , BD0T(1O,3) 

X06 ( 120 ) 



r*8 

Substitution symbol for 
DIN (10,3), DIND0T (10,3), 
D0GT ( 10, 3 ) , D0UTDT (10,3) 
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SUBROUTINE NAME: REA DC 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To prepare input words to be read by READIN 

using the subroutine SETUP 

CALLING SEQUENCE: Call READC 

LABELED COMMONS USED: IVARYP, EVARYP, BCNVRT, ISIMDT, IRHIST, 

XDBUG , IR0SIE, C/UFCN, CRjZfSIE, ICNC0R, TWRIT 

SUBROUTINES CALLED: SETUP 


GLOSSARY OF VARIABLES: All variables defined in labeled common 



SUBROUTINE NAME : READH 
LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To read flight data sets and store the 

pertinent data in a data array 

CALLING SEQUENCE: Call READH 

LABELED COMMONS USED: CSDATA, IRHIST, XDBUG 


SUBROUTINES CALLED: ACNVRT 

GLOSSARY OF VARIABLES: All variables defined in labeled common 
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SUBROUTINE NAME: r/sy 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: This subroutine evaluates the success or failure 

of a trail and provides stopping logic 

CALLING SEQUENCE: Call R/SY 

LABELED COMMONS USED: C/R2fSY, C/UFCN, CRgfelE, CR/SY, ICNC/R, 

IRQ S IE, PARRAY, RQSY1 

SUBROUTINES CALLED: GRAM 


GLOSSARY OF VARIABLES: All variables defined in labeled common 



SUBROUTINE NAME: SRINIT 


LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91*95 (OS) 

PURPOSE: To set up an indexing array to extract the appropriate 

simulator output which is to be compared with real data 
to determine performance 

CALLING SEQUENCE: Call SRINIT 

LABELED COMMONS USED: ISDARR, ISIMDT, XDBUG . IRjZfelE. IMAIN1, 

ANTENA 


SUBROUTINES CALLED: No subroutines 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

IS1(3) 



1*4 

Substitition symbol for 
IS ATT ( 3 ) 

IS2(3) 



1*4 

Substitution symbol for 
ISATTR( 3 ) 

IS3(2) 



1*4 

Substitution symbol for 
ISPHID, IS PHI 

IS4(120) 



1*4 

Substitution symbol for 
ISA (10,3), ISAD0T(1O,3), 
ISB( 10, 3 ) , ISBD/T(10,3) 

IS5(120) 



1*4 

Substitution symbol for 
ISDIN( 10, 3 ) , ISDIND(10,3), 
ISD/UT(10,30), ISD/UD(10,3] 
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SUBROUTINE NAME: UFCN (ICHECK, TIME) 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: To compute the performance criteria for each data 

channel and the total performance criteria for 
the system 

CALLING SEQUENCE: Call UFCN (ICHECK, TIME) 

LABE LED COMMONS USED: IRHIST, CSDATA, SIMDAT, XDBUG , IRjZfSIE, 

PARRAY, PR/S IE, C0UFCN, i/UFCN, RMAIN1 


SUBROUTINES CALLED: No subroutines 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 


MATH 

SYMBOL 


DIMENSION 


TYPE 


DESCRIPTION 


ICHECK 


1*4 


TIME 


sec 


R*8 


ICHECK 1 Initialize and 
preset values at start of 
trial 

ICHECK 2 Compute indi- 
vidual performance criteria 
for each data channel at 
each data point in a trial 

ICHECK 3 Compute total 
performance criteria at 
the end of a trial 

Time in seconds from the 
start of the year 
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SUBROUTINE NAME : WRITEP 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91.95 (OS) 

PURPOSE: To write standard corrector output 

at the end of each trial 

CALLING SEQUENCE: Call WRITEP 

LABELED COMMONS USED: C/UFCN, PR0SIE, CR/SY, IR/SIE, IRHIST, 

icnc/r 


SUBROUTINES CALLED: No subroutines 


GLOSSARY OF VARIABLES: All variables defined in labeled common 
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SUBROUTINE NAME: WRITER (ICS, IOS ) ENTRY READR 

LANGUAGE: FORTRAN IV 

MACHINE: IBM 360/75,91,95 (OS) 

PURPOSE: WRITER (ICS, Ijzfs ) Compute remaining CPU and 10 time 

and write corrector restart tape. READR to read corrector 
restart tape 

CALLING SEQUENCE: Call WRITER (ICS, i/s ) or Call READR 

LABELED COMMONS USED: C/r/SY, CSDATA, ROSYl, XDBUG, ISDARR, 

ICRPIN, PARRAY, CR/SY 

SUBROUTINES CALLED: REMTIM MAXO EXIT 


GLOSSARY OF VARIABLES: 


FORTRAN 

SYMBOL 

MATH 
. SYMBOL 

DIMENSION 

TYPE 

DESCRIPTION 

ICS 



1*4 

CPU time left at the 
start of the trial 

ijzfe 



1*4 

10 time left at the start 
of the trial 

ICST 



1*4 

CPU time left at the end 
of the trail 

IjZfeT 



1*4 

10 time left at the end 
of the trial 
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3 .2 Definition Of Labeled Common 

The new labeled common Is defined In the following paragraphs . 
The Interaction between new labeled common and subroutines Is 
shown In Figure 5. A description of the modifications to existing 
common is given below. 

C0M0N/PRC0^/ 

The array ST0RE (10, 30) was reduced to this size to 
conserve space. 

cf6m0N/c0yis0L/ 

The subscripted internal force coefficients and temperature 
coefficients were removed to conserve space. 

BLANK COMMON 

COlY0YrON/CO(7OOO) 

This common was labeled COMMON/C0RBIT/CO(7OOO) for ease of 
insertion in segment 4. 

C0MM0N/I0(300) 

This common was labeled COMMON/I i 0RBIT/I / 0(3OO) for ease of 
insertion in segment 4. 
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Figure 5 - Labeled Common - SUBROUTINE Interaction 
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N AME : C OMMON/BCNVRT/bSTIF 



-36 



NAME : COMMON/CSDATA/A ( 5000 ) 
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NAME: COMMON/ C/TJFCN/TWRDD , CWEGHT(20) J> CN0RM(2O) j) U, RESDUL(20) 
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NAME : COMMON/CjZfRjZfUT/PHSAVE ( 12 ) 




NAME: COmON/C0R0SY/C(12,12), D(12), DP (12), ULAST, USTAG, UPREV 
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EVABD(3) 0.0 DO feet R*8 Initial variation for element dis- 

placement initial conditions A (10, 3) 
: B(10,3), DIN ( 10, 3 ) j DjZfUT(l0,3), 

subscripted to mode 



NAME : COMMON/EVARYP/ (Cont ' d ) 
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differential TDIS 



NAME: C OMMON/EVARYP/ ( C on t ' d ) 
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NAME: COMMON/EVARYP/( Cone I'd) 
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bending stiffness coefficient BSTIF 



NAME: C0MM0N/ISIMDT/ISATT(3), ISATTR(3), ISPHID, IS PHI, ISA(10,3), ISAD0T( 10, 3) , ISB(10,3), 

ISBD0T (10,3), ISDIN (10,3), ISDIND(10, 3) , ISD0UT( 10, 3) , ISD0UD(1O,3) 
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NAME: C0MM0N/lCRPIN/lPAR(12,2) 
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NAME: COMMON/lSDARR/lSAR (20,2) 
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NAME: COMMON/lCNCjZfR/lC0RIC, NPARAM 
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NAME: C0MM0N/IRH1ST/NDW0RD, ITWRD(2) j IUW0RD(20), IDCHAN 
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placement initial conditions A (10, 3) 
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initial displacement DjZfUT ( 10, 3 ) 



NAME: COMMON/lVARYP/ (Cont ' d ) 
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sion spring constant, ZK1D 
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NAME: C OMMON/lR^S IE/MAXTRY , NWT, MMAG 
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NAME: COMMON/l0UFCN/lPU(2O) , IEUFCN 
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NAME : C OMMON/PARRA Y/P (12), E ( 12 ) 
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IREIN 0 regular start 
IREIN 1 restart 



NAME: COMMON/XDBUG/IDBUG, IFLITE, IZEF3, IZEF4 
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3 .3 Definition of Input 


The new and modified input symbols are defined in the following 
paragraphs . 

3.3.1 Generation of Test Plight Data Sets 


FORTRAN MATH 

SYMBOL SYMBOL DESCRIPTION UNITS 


NHISTS Control word for generating Standard Integer 

Test Plight Data Set. N.D. 

NHISTS 0 No data set generated 
NHISTS 1 Standard data set generated 
depending on IN0PT 

NHISTS 1 IN^PT 1 each data 
Record contains 5 double precision 
words 

Date time PSI1 PHIl THET1 


NHISTS 1 IN^PT 2 each data 
Record contains 14 double precision 
words 

DATE TIME ALPAE BETAE GAMAE 
A(l,l) A(2,l) A(3,l) A(4,l) B(l,l) 

B(2,l) B(3,l) B(4,l) PHILD 

IFLITE Control word for generating General Integer 

Test Plight Data Set N.D. 

IFLITE 0 No data set generated 
IFLITE 1 General data set generated 
by writing out complete 
normal output data record 
from GPRINT. Each data 
record contains store(300) 

Note: NHISTS 0 When IFLITE 1 

and IFLITE 0 When NHISTS 1 



3.3.2 General Control Words 


FORTRAN 

SYMBOL 

J0PTIM 
HP ARAM 


MATH 

SYMBOL DESCRIPTION UNITS 


Control word for corrector run 
J0PTIM 0 Normal simulator run 
jpPTIM 1 Corrector run 

Control word specifying the number 
of parameters to be varied in 
optimising the fit to flight data 


NPARAM 1 One Parameter 
NPARAM 2 Two Parameters 
Etc . 

Note: 1 ^ NPARAM ^ 12 

Control word specifying the type 
of performance criterion to be 
used 


J- I 
J = 2 
J = 3 


J = -f 


V‘li /V 

6 J- 

v-LL > s il 

i 6 




Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 
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FORTRAN 

SYMBOL 

MAXTRY 

NPRIN 


IDBUG 


MATH 

SYMBOL DESCRIPTION 


The maximum number of trials to 
be performed in optimising the 
fit to the flight data 


Control word to suppress standard 
simulator output for every trial 
For instance; NPRIN 10 gives 
regular simulator output at 
the first trial and every 
10th trial thereafter 

NPRIN 1 gives regular simulator 
output for every trial 


Control word for corrector 
DEBUG output . Not normally 
used 


IDBUG 0 No Output 
IDBUG 1 Debug, Output 


UNITS 


Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 


- 67 - 



3.3.3 Specification of parameters to be varied 

The user specifies the order in which parameters are to 
be varied by assigning the appropriate integer to the parameter 
control word; 1 is varied first; 2 is varied second; etc. Note: 
no more than 12 parameters can be varied in a single corrector 


run. 




FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DESCRIPTION 

UNITS 

WATT( 3 ) 


Control word for varying initial 
Conditions for attitude angles 

Integer 

N.D. 


IN0PT 1 WATT(l) - psil 
WATT (2) -THET1 
WATT(3) - PHI1 

IN^PT 2 WATT(1] - ALFAE 
IVATT (2) - BETAE 
WATT ( 3 ) - GAMAE 

WATTH(3) Control word for varying initial Integer 

conditions for attitude angular N.D. 

rates 


IN$PT 1 WATT R | 

1 1 > 

1 - 0MEGI 

(1) 

IVATTR 1 

2 

1 - 0MEGI 

2) 

WATTR 1 

[ 3 ! 

1 - j^MEGi 

[3) 

IN^PT 2 WATTR 1 

[1] 

1 - jZfMBCl 

fl) 

IVATTR 1 

2 

1 - 0mc\ 

2) 

WATTR 1 

3 

1 - 0MBCI 

(3) 


IVBELD Control word for varying damper Integer 

frame angle BETLD N.D. 
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FORTRAN MATH 

SYMBOL SYMBOL 

DESCRIPTION 

UNITS 

IVGALD 

Control word for varying the 
damper frame angle GAMLD 

Integer 

N.D. 

IVPHI 

Control word for varying the 
initial condition on the damper 
angular displacement PHILD 

Integer 

N.D. 

IVPHID 

Control word for varying the 
initial condition on the damper 
angular velocity DPHILD 

Integer 

N.D. 

iva(io,3) 

Control word for varying the 
initial conditions on the element 
A displacements 

Integer 

N.D. 

IVADp'r(10,3) 

Control word for varying the 
initial conditions on the element 
ADOT velocities 

Integer 

N.D. 

IVB(10,3) 

Control word for varying the 
initial conditions on the element 
B displacements 

Integer 

N.D. 

IVDIN( 10, 3 ) 

Control word for varying the 
initial conditions on the damper 
element DIN displacements 

Integer 

N.D. 

IVDIND( 10,3) 

Control word for varying the 
initial conditions on the damper 
element DIND0T velocities 

Integer 

N.D. 

IV D0UT ( 10, 3 ) 

Control word for varying the 
initial conditions on the damper 
element DjZfUT displacements 

Integer 

N.D. 

IVD0tJD(lO,3) 

Control word for varying the 
initial conditions on the damper 
element D^UTDT velocities 

Integer 

N.D. 

IVALPK(IO) 

Control word for varying the 
element frame EULER angle 
ALFAEK 

Integer 

N.D. 
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FORTRAN MATH 

SYMBOL SYMBOL 

DESCRIPTION 

UNITS 

IVBETK(IO) 

Control word for varying the 
element frame euler angle 
BETAEK 

Integer 

N.D. 

IVGAMK(IO) 

Control word for varying the 
element frame euler angle 
GAMAEK 

Integer 

N.D. 

IVBDj2fr(lO,3) 

Control word for varying the 
initial conditions on the element 
BDOT velocities 

Integer 

N.D. 

IVEM0B(1O) 

Control word for varying 
Young's modulus,, EM0DLS, for 
the element 

Integer 

N.D. 

IYTDIS(IO) 

Control word for varying 
the temperature differential 
across the element TDIS 

Integer 

N.D. 

IVZLO(IO) 

Control Word for varying 
the element length ZLO 

Integer 

N.D. 

IVSKjZfA ( 10, 3 ) 

Control word for varying the 
no load offset for the element 
in the A direction SK$A 

Integer 

N.D. 

IVSK0B( 10,3) 

Control word for varying the 
no load offset for the element 
in the B direction SKjZfe 

Integer 

N.D. 

IVDCAY 

Control word for varying the 
exponential decay constant in 
the magnetic hysteresis damper 
simulate. DECAY 

Integer 

N.D. 

IVZMD0 

Control word for varying the 
saturation torque in the magnetic 
hysteresis damper simulation 
ZMD0 

Integer 

N.D. 

IVZMDB 

Control word for varying the 
initial magnetic torque in the 
magnetic hysteresis damper simu- 
lation. ZMDBjZf 

Integer 

N.D. 

IVZK1D ' 

Control word for varying the 
damper suspension spring constant 
ZK1D 

Integer 

N.D. 
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FORTRAN MATH 


SYMBOL 

SYMBOL 

DESCRIPTION 

UNITS 

IVZK2D 


Control word for varying the 
damper suspension stop spring 
constant ZK2D 

Integer 

N.D. 

IVPHIS 


Control word for varying the 
damper suspension stop angle 
PHIS 

Integer 

N.D. 

IVCNV 


Control word for varying the 
damper suspension viscous torque 
coefficient CNV 

Integer 

N.D. 

IVDT0# 


Control word for varying the 
temperature differential across 
all elements simultaneously DT00 

Integer 

N.D. 

IVSTIF 


Control word for varying the 
bending stiffness of all elements 
simultaneously BSTIF 

Integer 

N.D. 

IVCDMP(3,10) 


Control word for varying the 
damping factor for antenna motions 
CDAMP 

Integer 

N.D. 

3.3.^ Initial Variations for Parameters 


The user 

must 

specify an initial variation step for 

the 

parameters to 

be varied in a corrector run. 


FORTRAN 

SYMBOL 

MATH 

SYMBOL 

DESCRIPTION 

UNITS 

EVATT 


Initial step size for variations 
in attitude angle initial condi- 
tions . Must be specified if 
IVATT(l-3) non zero 

Degrees 

EVATTR 


Initial step size for variations 
in attitude angular rate initial 
conditions . Must be specified 
if IVATTR (1-3) non zero 

Degrees/si 

EVPHIL 


Initial step size for variations 
in damper reference frame Euler 
angles BETLD GAMLD PHILD. Must 
be specified if IVBELD, IVGALD, 
or IVPHI non zero 

Degrees 
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FORTRAN 

SYMBOL 

EVPHID 
EVABD ( 3 ) 


EVABDD(3) 


EVALPK 

EVBETK 

EVGAMK 


MATH 

SYMBOL DESCRIPTION 

Initial step size for variation 
in damper motion angular rate 
initial condition DPHILD. Must 
be specified if IVPHID non zero 

Initial step size for variations 
in element displacement initial 
conditions. Must be specified 
if IVA (1-10,3), IVB( 1-10,3), 

IVDIN( 1-10, 3 ) , or IVD0UT(1-1O,3) 
are non zero 

EVABD ( 1 ) for A(l-10,l), B(l-10,l), 

DIN( 1-10,1), D0UT( 1-10, 1 
EVABD ( 2 ) for A(l-10,2), B(l-10,2), 

DIN( 1-10,2), D0UT(1-1O,2 
EVABD ( 3 ) for A(l-10,3), B(l-10,3), 

DIN ( 1-10 , 3 ) , DjZfUT ( 1-10, 3 

Initial step size for variations in 
element velocity initial conditions . 
Must be specified if IVAD0T( 1-10, 3 ) , 
IVBD0T (1-10,3 ) , IVDIND( 1-10, 3 ) , or 
IVD0UD(1-1O,3) are non zero 

EVABDD ( 1 ) for AD0T( 1--10, 1 ) , BD0T(1-1O 
DIND0T(1-1O,1), 

D0UTDT( 1-10, 1 ) 

EVABDD ( 2 ) for AD0T( 1-10, 2 ) , BD0T(1-1O 
DIND0T(l-lO y 2), 

D0UTDT( 1-10,2 ) 

EVABDD (3 ) for AD0T(1-1O,3), BD0T(1-1O 
DIND0T(1-1O,3), 

D0UTDT ( 1-10, 3 ) 

Initial step size for variation in 
element reference frame Euler 
angle ALFAEK. Must be specified 
if IVALPK (1-10) non zero 

Initial step size for variation in 
element reference frame Euler 
angle BETAEK. Must be specified 
if IVBETK(l-lO) non zero 

Initial step size for variation in 
element reference frame Euler 
angle GAMAEK. Must be specified 
if IVGAMK(l-lO) non zero 


UNITS 

Degrees/sec 

Feet 


Feet/sec 

1 ) 

2 ) 

3) 

Degrees 

Degrees 

Degrees 
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FORTRAN 

SYMBOL 

EVEM/D 

EVTDIS 

EVZLO 

EVSK0(3) 

EVDCAY 

EVZMD/ 

EVZMDB 

EVZK1D 


MATH 

SYMBOL DESCRIPTION 

Initial step size for variation 
in Young's modulus of element 
EM/DLS . Must be specified if 
rVEMjZfo(l-lO) non zero 

Initial step size for variation 
in the temperature differential 
across the element TDIS . Must 
be specified if IVTDIS ( 1-10) 
non zero 

Initial step size for variation 
of element length ZLO . Must be 
specified if IVZLO(l-lO) non zero 

Initial step size for variation 
of element no load offsets SK/A, 

SKjZfe. Must be specified if 
IVSKjZfA ( 1-10, 3 ) or IVSK/B( 1-10, 3 ) 
are non zero 

EVSK/(l) for SK0A( 1-10,1), SK#B(1-10, 
EVSKjZf ( 2 ) for SK0A (1-10,2 ) , SK#B(1-10, 
EVSEYS) for SK/a( 1-10, 3) , SK#B(1-10, 

Initial step size for variation 
of the exponential decay constant 
in the magnetic hysteresis damper 
simulation. Must be specified if 
IVDCAY non zero 

Initial step size for variation 
of the saturation torque in the 
magnetic hysteresis damper simu- 
lation. Must be specified if 
IVZMD/ non zero 

Initial step size for variation 
of the initial magnetic torque 
in the magnetic hysteresis damper 
simulation. Must be specified if 
IVZMDB non zero 

Initial step size for variation 
of the damper suspension spring 
constant ZK1D. Must be specified 
if IVZK1D non zero 


UNITS 

2 

Pounds/inch 


N.D. 


Feet 


Feet 


1 ) 

2 ) 

3) 

N.D. 


Foot pounds 


Foot pounds 


Foot pounds/ 
rad 
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FORTRAN MATH 

SYMBOL SYMBOL 

EVZK2D 


EVPHIS 


EVCNV 


EVDT// 


EVSTIF 


EVCDMP ( 3 ) 


DESCRIPTION 

Initial step size for variation 
of the damper suspension stop 
spring constant ZK2D. Must be 
specified if IVZK2D non zero 

Initial step size for variation 
of the damper suspension stop 
angle PHIS . Must be specified 
if IVPHIS non zero 

Initial step size for variation 
of the damper suspension viscous 
torque coefficient CNV. Must 
be specified if IVCNV non zero 

Initial step size for variation 
of the temperature differential 
across all elements simul- 
taneously. Must be specified 
if 1VDT00 non zero 

Initial step size for variation 
of the bending stiffness for 
all elements simultaneously. 

Must be specified if IVSTIF 
non zero 


UNITS 

Foot-pounds/ 

rad 


Degrees 


Foot-pounds/ 

sec 


o 


F 


N.D. 


Initial step size for variation N.D. 

of element damping coefficient 

EVCDMP(l) for CDAMP(l,l-10) 

EVCDMP (2 ) for CDAMPf 2,1-10] 

EVCDMP(3) for CDAMP( 3, 1-10) 

Must be specified if IVCDMP(3* 1-10) 
non zero 


3*3.5 Reading of Flight Data Sets 

The input described below is required in the reading of 


general flight data sets only . Standard flight data sets, defined 
in the format by NHISTS and IN/PT, will be read automatically by 
the program without special input . 

The general flight data set is assumed to be a series of 
sequential data records with no header record. Each data record 
is assumed to be a vector of Real*8 words. Each data record is 
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assumed to contain the date (YYMMDD.) year, month, day, and 
time of day (HHMMSS . ) hour, minute, second. 

FORTRAN MATH 

SYMBOL SYMBOL DESCRIPTION 

ISTACK Control word for presetting data 

for the reading of standard 
flight data sets 
ISTACK 0 standard flight data 
ISTACK non zero general flight 
data set 


The following input words are required only for ISTACK 


non zero: 
NDWORD 


ITWRD( 2 ) 


IDCHAN 


IUW0RD(20) 


The number of Real*8 words in 
each data record 

For general flight date sets 
generated for testing purposes 
with the program NDWORD 300 

The location in the data record 
of the two time words 

DATE = YYMMDD . 

HMS = HHMMSS. 


ITWRD (1 
ITWRD ( 2 


is location of DATE 
is location of HMS 


For general flight data sets gen- 
erated for testing purposes with 
the program 

ITWRDfl) 1 
ITWRD (2 ) 2 

The number of data channels to be 
used in the evaluation of the 
performance criteria. The maximum 
number of data channels is 20 

The location in the data record 
of the various data channels for 
instance, with a general flight 
data set generated for testing 
purposes with the program and 
INjZfPT 2 

ALFAE in data channel 
BETAE in data channel 
GAMAE in data channel 


IUWjZfRDfl 
IUW0RD 2 
IUW0RD ( 3 



the 


UNITS 

Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 
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FORTRAN 

MATH 



SYMBOL 

SYMBOL 

DESCRIPTION 

UNITS 

TWRDD 


Test word for recognizing missing 
data 

N.D. 


3.3.6 Simulator Data Channels for Performance Criteria 

The input required to assign the appropriate data channel 
to simulator output in order to construct the performance criteria 
is described below. In the case where the flight data set is 
of a standard format consistent with IN0PT the simulator output 
is assigned automatically by use of ISTACK 0. For standard 
flight data sets the data channel assignments are as follows: 


CHANNEL NO. IN0PT 1 INOPT 2 


1 


PSI1 

ALFAE 

2 


PHI1 

BETAE 

3 


THET1 

GAMAE 

4 


UNUSED 

A(l,l) 

5 


1! 

A 2,1 

6 


M 

A (3,1) 

7 


n 

A (4,1) 

8 


M 

b(i,i) 

9 


II 

B 2,1 

10 



B 3,1 

n 


It 

b(4,1) 

12 • 


M 

PHILD 


No additional input is required to assign the above formats using 
standard flight data sets other than ISTACK 0 and IN0PT 1 or 2 . 

In the use of general flight data sets it is necessary to 
provide input to assign the appropriate simulator output quantities 
to the proper data channel. For general flight data sets the 
control word ISTACK must be non zero . The assignment of the 
flight data to particular data channels is made by the input of 
the vector IUW0RD(20) in extracting the data from the general 
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flight data set . The assignment of the simulator output data to 
construct the performance criteria must be compatible with the 
data channels assigned for the flight data. 


FORTRAN 

SYMBOL 

ISATT(3) 


NOTE: J = 

isattr( 3 ) 

IS PHI 
ISPHID 
ISA (10,3) 
ISAD0T(1O,: 

ISB(10,3) 


MATH 

SYMBOL DESCRIPTION 


Control word to assign simulator 
output to appropriate data channel. 

The meaning of ISATT(3) is modi- 
fied by INJ3PT 

IN0PT 1 

ISATT(l) J assigns PSI1 to channel J 
ISATT(2) J assigns THET1 to channel J 
ISATT( 3 ) J assigns PHI1 to channel J 


IN0PT 2 

ISATTfl] J assigns ALFAE to channel J 
IS ATT (2) J assigns BETAE to channel J 
IS ATT ( 3) J assigns GAMAE to channel J 


0 data not assigned 


Control word to assign simulator 
output to appropriate data channel. 

The meaning of ISATTR(3) is modi- 
fied by IN0PT 

INjZfPT 1 

ISATTR(I) J assigns 0MEG ( I ) to channel 
IN0PT 2 

ISATTR(l) J assigns 0MBC(l) to channel 


Control word to assign simulator 
output damper angle PHILD to appro- 
priate data channel 

Control word to assign simulator 
output damper angular rate DPHILD 
to appropriate data channel 

Control word to assign simulator 
output element displacement A (10,3) 
to appropriate data channel 

Control word to assign simulator 
output element velocity AD^T(10,3) 
to appropriate data channel 

Control word to assign simulator 
output element displacement B(10, 3) 
to appropriate data channel 


UNITS 

Integer 

N.D. 


Integer 

N.D. 


J 


J 

Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 


Integer 

N.D. 
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FORTRAN 

SYMBOL 


MATH 

SYMBOL 


DESCRIPTION 


UNITS 


ISBD0T(1O,3) 

ISDIN(10,3) 

ISDIND( IOj 3 ) 

ISD0UT( 10, 3 ) 

ISD0UD( 10, 3 ) 


Control word to assign simulator Integer 

output element velocity BDjZfT(10,3) N.D. 

to appropriate data channel 

Control word to assign simulator Integer 

output damper element displacement N.D. 

DIN(10,3) to appropriate data 
channel 

Control word to assign simulator Integer 

output damper element velocity N.D. 

DINDjZfT(10,3) to appropriate 
data channel 

Control word to assign simulator Integer 

output damper element displacement N.D. 

D0UT( 10,3) to appropriate data 
channel 

Control word to assign simulator Integer 


output damper element velocity Dj2fUTDT (10, 3) 
to appropriate data channel 


3.3.7 Performance Criteria 

The performance criteria is constructed using the residuals 
from each data channel either directly or weighted and normalized. 
The input symbols used to provide weighting and normalization are 
described below: 


FORTRAN 

MATH 



SYMBOL 

SYMBOL 

DESCRIPTION 

UNITS 

CWEGHT(20) 


Weighting factor for residuals 
subscripted to the data channel. 
These are used only with 
J = 3 or 4 

N.D. 

CN0RM(20) 


Normalization constant for 
residuals subscripted to the 
data channel . These are used 
only with J = 3 or 4 

Units of 
Data 
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3.3.8 Input for Parameter Search Operation 

In the process of minimizing the performance criteria 
certain constants are required in computing parameter variations 
and in testing the performance criteria. The input symbols for 
these constants are described below. However, the use of preset 


values is 
FORTRAN 

recommended 

MATH 

for routine operation. 


SYMBOL 

SYMBOL 

DESCRIPTION 

UNITS 

C ALPHA 

oC 

Factor in selecting next step 
size of parameter after a success- 
ful trial 

OC > 1 

N.D. 

CBETTA 

/3 

Factor in selecting next step 
size of parameter after an un- 
successful trial 

/3 < i 

N.D. 

CGAMMA 

r 

Factor in selecting next step 
size of parameter after a success- 
ful stage 

N.D. 

RATU 

R 

Stopping criteria constant ratio 
of difference in three successive 
stage performance criteria 

0.0 < R < 1.0 

N.D. 

DEL 

D 

Band in first completion criteria 
0 < D < 1.0 

N.D. 

T0L 


Tolerance for determining success 
of a given trial 

N.D. 



3 . 3.9 Restart Procedure 

Normal corrector operation to determine best fit parameters 
will require relatively large amounts of machine time. The 
program is designed to check remaining CPU and 10 time at the 
end of each trial. When either remaining time approaches the 
maximum time for any trial a restart tape is written so that the 
search can be continued from that point automatically. 

FORTRAN MATH 

SYMBOL SYMBOL DESCRIPTION UNITS 

Control word to read a restart Integer 

tape to continue a corrector run. N.D. 

IREIN 0 not a restart 

IREIN 1 read restart tape 

NOTE: To use a restart tape the 

original data deck must 
be used with the addition 
of IREIN 1 

3 .4 Description of Output 

The normal corrector output is printed at the end of each 

trial. The first line consists of three control words which 

describe the state of the search procedure . 

NSTAGE The number of completed searches through 
the parameter space. 

NTRIAL The number of the trial. This counter is 
incremented after each trial. 

NSUCCESS The number of consecutive trials in any 
parameter direction in which the per- 
formance criteria was reduced. If 
NSUCCESS is non zero the current trial 
reduced the performance criteria. 

The vector of parameter values for the current trial is printed 

out under the heading parameter values. The total performance 

criteria for the current trial is printed and the individual 

residuals of each data channel are recorded. 
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3 .5 Overlay Structure 

This section describes the organization of the program and 
the overlay structure. The program consists of four modules. 

A description of each module by size and composition is presented 
and the overlay tree structure is shown in Figure 6. A job 
control language listing is also included. 
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1. Root Segment A, Size 


197 K 


HAIY^UT 

IDATE1 

RP^L5 

XIN2 

GLITCH 

pi^tit 

ARTLU1 

IN0UT- 

NUM 

TYPSET 

SUNANG 

EVARYP 

CR0SY 

ANTENA 

XDBUG 

MOMENT 

PL0T 

2 . Segment 

READER 

READH 

3 . Segment 

ADCAL 

S^LAR 

PASTR 

NXNS^L 

^MCL 

VISCUS 

WHEELS 

G PRINT 

UFCN 

RP^L9 

STAPES 

4 . Segment 

ADM4RK 

Output 

READE 

C0RBIT 


SATP^S 

ZSPINR 

RP00L3 

SWITCH 

CSTVAL 

RVISCS 

CGMP62 

OTIME 

LAOOOOOO 

ADCN 

ICRPIN 

CI^SIE 

PRC0M 

IFEMRS 

SIMDAT 

CS0LAR 

AICRT3 


ITEMP2 

XIN4 

RP$0L1 

RTDIST 

CjZfoPLY 

RMGNTC 

DR0TJT 

STD PND 

MIN 

BCDWD 

PARRAY 

BCNVRT 

ENERGY 

CEFMRS 

RFIELD 

C0NSTS 

ENERGY 


0UT0NE RMAIN2 
JBTEST ITEMP 
RDAMPR IjZfDPLY 
IFASTR IMAIN2 
CDMPER ADREAL 
ACNVRT ARFDT2 
HI HMSGHJT 

TCNVRT YR]V(0DA 
ANUMARG DECK1 
CSDATA PR0SIE 
IR^SIE ISIMDT 
TWRIT C^UFCN 

C0MBCR R0SY1 
RP00L2 IP00L1 
RATTDE PWHEEL 
CFNALP CCNVRT 
PLTND 


DATGflJT 

CCMBNZ 

HEAD22 

XIN1 

DEPL0Y 

N^DER 

CHECK 

multm 

FVAL 

SHFT2V 

LltfUFCN 

IVARYP 

ISDARR 

RMAIN1 

C0R0SY 

ECNSTS 

PL0TF 

B, Size 

READZ 


C, Size 

ADMIMP 

C0MBNZ 

FIELD 

MYPMAT 

PINTRP 

TPREAD 

XACM 

G0UT 

WRITER 

^UTTW^ 

RP^L4 

D, Size 

ASINR 

^RBEL 

SUB3 

I^RBIT 


23 K 
REA DC 


230 K 

AIRDRG 

DAMPER 

FINDGB 

MATMPY 

PLH 

SUNLIN 

XYZPLH 

Rj2fSY 

S,0UT 

ELMNTS 

RP00L6 

106 K 

ATANQR 

N^INTP 

SUB4 

I^UTPT 


BL0CK4 


ARANGE 

DARC0S 

FINDS 

MAGNTS 

REL0C 

SUNDEP 

ZSPIN 

C0NVRT 

CANTNA 

VARBLS 

^UTTHR 


C^MPEX 

MAIN2 

SUB5 

AC0SR 


EISUBK 

XIN3 

RFASTR 

RSUNCL 

CFINDX 

ADDAM 

EL0RJT 

SIVRBL 

SBLOGOR 

C^RE 

ICNC/R 

IRHIST 
IMA INI 
CRATI0 
CFNDGB 
IDFRMV 


CRINIT 


ATTUDE 

DARSIN 

FNDALP 

IS PIN 

REfyZfC2 

SUBC^N 

GRAM 

LSARAY 

DERFND 

RP$0L8 


C^NVIN 
JULDAY 
TRANS P 


SRINIT 


B0UNDS 
DEREQl 
FNDGAM 
INVERT 
RTS IDE 
S INPUT 
GINPUT 
LDARAY 
C0MALP 
C0MSjZfL 


DEREQ 

GETTAP 

ZER^ 


ARFIC 


CKINPT 

DSAVE 

FNYBCM 

HAG 

SET 

SETVAL 

XFIND 

WRITEP 

RP00L7 

^UTF0R 


FMERIS 

FINDX 

AFMRIS 
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JOB CONTROL LANGUAGE LISTING 


Column 

1 


//YVTLBOOl JOB ( P20101861F, P, NA0103, 004005 ) , 077, MSGLEVEL= (1,1) 
//LINK EXEC LINKGO, PARM. LINK= * OVLY, LET, LIST, XREF, MAP ' , 

// REGION. G0=480K,C0ND=(382,EQ) 

//LINK. TAPE LIB DD DSN=PLB, V0Ip=SER=TD6052, LABEL= ( 1, BLP ) , 

// DISP= ( OLD, DELETE ) 

//LINK. SYS LIN DD * 


* 

* 


PUT AVCOPAC OBJECT DECK HERE 


INCLUDE TAPELIB 

INSERT LA000000,MAIN,ADDAM,ACNVRT,ARFDT2,ARTLU1, CHECK, C0MP62 
INSERT DROUT,ELOUT, FVAL, HI, HMSOUT, INOUT,MULTM, PLOT, PLOTF, QTIME 
INSERT STDPND, STVRBL,TCNVRT,YRMODA, AICRT3, IDFRMV, PLTND 
INSERT ADREAL, ANTENA, CCMBNZ, CCNVRT, CDMPER, CEFMRS, CFINDX 
INSERT CFNDGB, CFNALP, CODPLY, CONSTS, CRATIO, CSOLAR, CSTVAL, DATOUT 
INSERT DEPLOY, ECNSTS,EISUBK, GLITCH, HAMOUT, HEAD22, ID ATE 1 
INSERT IFASTR, IFEMRS, IMAIN1, IMAIN2, IODPLY, IP00L1, ITEMP, ITEMP2 
INSERT JBTEST, MOMENT, NODER, OUTONE, PLOTIT, PRCOM, PWHEEL 
INSERT RATTLE, RDAMPR, RFASTR, RFIELD, RMAIN1, RMAIN2, RMGNTC 
INSERT RP00L1, RP00L2, RP00L3, RP00L5, RSUNCL, RTDIST, RVISCS 
INSERT SATPOS, SUNANG, SWITCH, XIN1,XIN2,XIN3, XIN4, ZSPINR 
INSERT NUM 
OVERLAY A1 

INSERT CRINIT, READH, SRINIT, READER, READZ, READC, BL0CK4, ARFIC 
OVERLAY A1 

INSERT ROSY, WRITE P, WRITER, GRAM, LDARAY, UFCN, LSARAY 
INSERT ADCAL, ADMIMP, AIRDRG, ARANGE , ATTUDE , BOUNDS, CKINPT 
INSERT COMBNZ, CONVRT, DAMPER, DARCOS,DARSIN, DEREQ1, DSAVE , FASTR 
INSERT FIELD, FINDGB, FINDS, FNDALP, FNDGAM, FNYBCM, GINPUT, GOUT 
INSERT GPRINT,HAG, INVERT, IS PIN, MAGNTS,MATMFY,MPYMAT 
INSERT NXNSOL, OMCL, PINTRP, PLH, RELOC, REL0C2, RTSIDE, SET 
INSERT SETVAL, SINPUT, SOLAR, SOUT, SUBCON, SUNDEP, SUNLIN, T PRE AD 
INSERT VISCUS, WHEELS, XA CM, XFIND,XYZPLH, ZSPIN 
INSERT COMALP, COMSOL, CANTNA, DERFND,ELMNTS, OUTTWO, OUTFOR 
INSERT OUTTHR, RP00L4, RP00L6, RP00L7, RP00L8, RP00L9 
INSERT ST A P OS , VA RBLS 
OVERLAY A1 


INSERT AFMRIS, ACOSR, ADM4RK, ASINR, ATANQR, COMPEX, CONVTN, DERE Q 
INSERT FMERIS, FINDX, GETTAP, JULDAY,MAIN2,N0INTP, ORBEL, OUTPUT 
INSERT READE, SUB3, SUB4, SUB5,TRANSP, ZERO, CORBIT, IORBIT, IOUTPT 
ENTRY MAIN 
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Column 

1 


JOB CONTROL LANGUAGE LISTING (Continued) 


GO. 


GO. 

GO. 


FTlOFOOl DD UNIT=(9TRACK, , DEFER), LABEL=(l,BLP), 

V0Ip=SER=TD5244 . DCB= (DEN=2, RECFM=VS, LRECB=7460, BLKSIZE=7464) 
DISP=( OLD, KEEP) 

TD5244 IS JPL EPHEMERIS TAPE 
FTllFOOl DD UNIT= (2400-4,, DEFER), LABEL=(l,BLP) , 

VOL=SER=POOL, DCB= (RECFM=VBS, LRECI^124,BLKSIZE=2484) , 

DISP= (OLD, PASS) 

FT12F001 DD UNIT=23l4, DISP= (NEW, DELETE ) , SPACE= (CYL, (2,2) ), 
DCB= ( RE CFM=VBS, LRE CL=4l2, BLKSIZE=4l6 ) 

FT15F001 DD LABEL=(1,BLP),UNIT=(9TRACK,, DEFER), 

DISP= (NEW, KEEP) ,V0L=SER=TD6212, 

DCB= (RECFM=VS, LRECL=796,BLKSIZE=800) 

TD6212 IS A FLIGHT DATA SET 

FT21F001 DD DSN=RESTART,LABEB= (l,BLP) ,UNIT= ( 9 T RACK, , DEFER) , 
V0L=SER=TD6401.DCB= (recfm=vs,lrecl=796,blksize=7964) , 

DISP=( NEW, KEEP) 

TD6401 is a restart tape 

SYSABEND DD SYSOUT=A 
DATA 5 DD * 
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3 .6 User's instructions 

The following pragraphs are intended to aid the user in the 
operation of the corrector option of this program. Whether the 
program is used in a simulation mode or a corrector mode, a 
complete set of data to define the satellite, orbit, initial con- 
ditions, and environment is required. All the control words for 
the corrector mode are preset so that corrector mode operation 
is not invoked. This means that a normal simulation mode run 
requires only the usual input . 

3.6.1 Writing Plight Data Sets 

Test flight data sets can be written with the program if 
desired. When writing a test flight data set, the program is 
operating in the simulation mode. The input control words 
NHISTS 1 and IFLITE 0 create standard flight data sets . The input 
control words WHISTS 0 and IFLITE 1 create general flight data 
sets . The input control words NHISTS and IFLITE are preset to 
zero. Since the corrector control words are preset to not invoke 
corrector operation, none of the corrector control words should 
be input to write a test flight data set. 

3.6.2 Corrector Operation 

There are three basic requirements for using the corrector. 

The flight data and simulator data must be compatible; the type 
of error function or performance criterion must be selected; and 
the parameters to be varied must be defined. A discussion of 
these aspects of corrector use is presented below. 

Flight data is read from a data tape and stored in common in 
a particular form for use by the corrector. The reading of flight 
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data is accomplished in the subroutine READH. The current program 
can read tapes created in either of the standard forms or tapes 


created in the general form described in Section 3.3. However, 
it is likely that the reading of real data tapes will require a 
modification to READH. In making this modification the structure 


of the resulting data storage in common should be strictly main- 
tained. The common storage of data has the following format. 

The array is in double precision. The first NDW0RD locations 
are used for the temporary storage of each data record as it is 
read from the tape. The subsequent locations in the array are 
used for sequential storage of the pertinent data from the data 
records read from the flight data tape. For each data record, 
the time in seconds from the start of the year is stored first. 


The data for the appropriate data channels is stored in the next 
IDCHAN locations . This format is repeated for each time point 


in the data span. The resulting common array is as follows: 


A(l) 

« > 

: i 

A ( NDW0RD ) j 

TSTART 
DATA ( 1 ) 

DATA (2) 

DATA ( 3 ) 

DATA (IDCHAN) j 

TIME 
DATA ( 1 ) 

DATA (2 ) 

DATA (3) 

DATA (IDCHAN) J 

TST0P 
DATA ( 1 ) 

DATA (2) 

DATA ( 3 ) 

DATA (IDCHAN) J 



Temporary storage 
for data records 
from flight data tape 


Stored data at 
time TSTART 


Stored data at 
intermediate time 


Stored data at 
time TST0P 
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In evaluating the performance criteria for the corrector* 
simulator data must be put into an array at each print time for 
the simulator. The simulator data must be put into the data 
channels in the same order as has been chosen for the flight 
data. This can be accomplished by the input described in 
Section 3.3. In the event that real flight data is not compatible 
with the normal data produced by the simulator* suitable adjust- 
ments in the program will be required. 

The selection of a particular type of performance criterion 
has no clear cut guidelines and only experience with each particu- 
lar situation can determine a "best" criterion. It is probably 
advisable to put all the data on some nondimens ional basis. For 
instance* if a criterion is being made up of both attitude angles 
and element displacements* the element displacements should 
possibly be nondimensionalized by weighting with the reciprocal 
of the antenna length. The best procedure may be to try several 
different criterion to see if any improvement in optimization 
can be observed. 

The mechanics of defining the parameters which are to be 
varied consists of specifying the appropriate parameter control 
words and initial variations as described in Section 3-3. The 
more difficult problem is to decide which parameters should be 
varied. The first step in fitting of a given data span is to 
optimize the initial state vector particularly with respect to 
rates. Once the initial conditions have been resolved* the inves- 
tigation of physical parameters of the system can be undertaken. 
There are no clear rules for selection of parameters. Selection 
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is more basically related to an evaluation of what parameters 
can explain the observed performance and an estimate of the 
inherent uncertainty in the knowledge of a particular parameter. 
The careful combination of these two concepts can lead to an 
estimate of a "most likely" set of parameters. 

3.6.3 General Instructions and Recommendations 

The following comments are intended as guidelines in using 
the corrector option. 

The internal orbit option should be used for all 
corrector operation. 

Intermittent printing of simulator output for 
corrector trials should be used to conserve output 
(NPRIN 10 is reasonable). 

Plotting of simulator output for corrector trials 
should not be used (iPLjZfT 0). 

When using a restart tape* the complete data deck 
for the original run must be used. 

When using a restart tape, the flight data tape is 
not required. 
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4.0 DETERMINATION OP ANTENNA POSITIONS PROM ATTITUDE DATA 
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An investigation was undertaken to evaluate the feasibility 
of determining definitive antenna positions on the basis of 
attitude data for a spacecraft of the RAE type. This study con- 
sisted of two distinct efforts . The first was an examination of 
a simplified linear formulation of the equations of motion. The 
second was the use of the corrector to demonstrate the problems 
of determining definitive antenna positions . These two areas 
are discussed in the following paragraphs . 

4.1 Simplified Equations of Motion 

The equations of motion for the RAE satellite are quite 
complex. For illustrative purposes it is desirable to consider 
a simpler set of equations of motion associated with an X boom 
configuration. The equations of motion are taken from Reference 4. 
The planar geometry of the situation is shown in Figure J. 



A 


i 


Figure 7 - Geometry of X Configuration 
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The equations of motion for pitch and in plane antenna displace- 
ments can be written in the following simplified form 

pitch equation of motion (linearized) 

I - M* f, + C a % + K v f, - 3 nl Lx, -rJ/? = m 2 

element in plane equations of motion (linearized) 

•» ** 

+ M, +• f, = Q, 

ft + c x*.fi + 4 fz = Q.* 

^ S'* "^3 ^3 ~ Q- s 

"V A+<fr 4 = < 2-4 

The element displacements associated with the generalized 
coordinates q^, q^, q^, and q^ are shown in Figure 8. 

An examination of these simplified equations of motion 
provides the essential information for evaluating the feasibility 
of determining definitive antenna positions from attitude data. 
The motions associated with the generalized coordinates q^ and qg 
are coupled to pitch but the motions associated with q^ and q^ 
have no effect on attitude in the linear approximation. Further- 
more, the motion associated with qg is only coupled to attitude 
by the fact that the antenna is bent by gravity gradient forces . 
From these observations, it is clear that the feasibility of 
determining definitive antenna position from attitude data is 
questionable. 
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fys — Ai + A2" A3 ~ A4. ^4 “ A 1 - Az“*~ A3 A4 

Figure 8 - Deformations for Generalized Coordinates 

4 .2 Illustrative Corrector Runs 

In order to demonstrate the problems associated with 
determining definitive antenna positions from attitude data, 
two corrector runs were made. In these runs, the initial condi- 
tions on the in plane antenna displacements were treated as 
parameters . The flight data for the illustrative cases was 
created by the program. This means that, as far as the corrector 
operation was concerned, the data was perfect and the simulator 
model was perfect. Therefore, the results obtained represent 
the best estimates of antenna positions that can be expected for 
the limited number of runs performed. 
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The true flight data initial conditions on in plane displace- 
ments for the two cases were the following: 

A 2 = 140.53 A 2 = -163.58' A 3 = - 162.92 A 4 = 111.55 

or q x = -74.42 q 2 = 578.58 q 3 = 27.98 q 4 = 29.64 

The initial estimates of the in plane displacements for the two 
cases were the following: 

A x = 140 .0 A 2 = -140.0 A 3 = -140.0 A^ = 140.0 

or q^ = 0.0 q 2 = 56 O.O q 3 = 0.0 q^ = 0.0 

The estimates roughly correspond to gravity gradient equilibrium 
displacements . 

The first corrector run considered a data span of sixteen 
minutes with points at thirty second intervals resulting in thirty- 
three data points for each data channel. The results of success- 
ful trials are given below for the first two stages . An unweighted 
squared error function was used. 


Trial 

No. 

137.5 

Aq. 

-140.0 

*4 

140 . 0 

Performance 

Criterion 

3 


-140.0 

63.055 

4 


132.5 

-140.0 

-140.0 

140.0 

54.688 

5 


122.5 

-140.0 

-140.0 

140.0 

39.774 

6 


102.5 

-140.0 

-140.0 

140.0 

17.204 

7 


62.5 

-140.0 

-140.0 

i4o.o 

0.996 

9 


62.5 

-135.0 

-140.0 

■i4o.o 

0.951 

12 


62.5 

-135.0 

-142,4 

140.0 

0.836 

17 


62.5 

-135.0 

-142.4 

139.4 

0.835 

End of 

First 

Stage q 1 = 

= -75.5, 

q 2 - 479.3, 

q 3 = • 

-69.5, q 4 = -84.3 

23 


64.9 

-135.0 

-142.4 

139.4 

0.829 

26 


64.8 

-136.3 

-141.8 

139.5 

0.833 

29 


64.8 

-136.6 

-142.4 

139-4 

0.809 

31 


64.8 

-136.6 

-142.5 

139.7 

0.813 

End of 

Second Stage q^ 

= -74.8, 

q 2 = 483.2 

, q 3 = 

-68.8, q 4 = -84. 
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The standard deviation of attitude error for the best performance 





0.809 

3TW 


0.090 


The second corrector run considered a data span of sixty 
minutes with points at thirty second intervals resulting in one 
hundred and twenty-one data points for each data channel. The 
results of successful trials are given below for the first two 
stages. An unweighted squared error function was used. 


Performance 


Trial No . 

At 


hr 

h 

Criterion 




3 

137.5 

-140.0 

-140.0 

140.0 

169.038 

4 

132.5 

-140.0 

-140.0 

140.0 

146.388 

5 

122.5 

-140.0 

-140.0 

140.0 

106.312 

6 

102.5 

-140.0 

-140.0 

140.0 

47.631 

7 

62.5 

-140.0 

-140.0 

140.0 

11.225 

11 

62.5 

-138.8 

-140.0 

140.0 

11.227 

14 

62.5 

-138.8 

-142.5 

i4o.o 

11.215 

16 

62.5 

-138.8 

-142.5 

145.0 

11.192 


End of 

First Stage q 1 : 

- -73.8, 

q 2 = 488.8, 

q 3 = ' 

■78.8, 

q 4 = -86.2 

22 

64.9 

-138.8 

-142.4 

144.8 


11.079 

24 

65.O 

-138.7 

-142.7 

145.4 


11.180 

27 

64.8 

-139.3 

-142.8 

145.5 


10.995 

28 

64.4 

-140.4 

-142.9 

145.8 


10.715 

29 

63.6 

-142.7 

-143.1 

146.3 


10.489 

32 

63.5 

-142.9 

-143.2 

145.1 


10.506 

End of 

Second Stage q^ 

= -75.9, 

q 2 = 495.7, 

' q 3 = 

-82.3, 

q 4 = -83.1 


The standard deviation of attitude error for the best performance 


is given by: 

<r“ 


10.489 

3(l2l) 


0.17 


The results of these two corrector cases illustrates the 
difficulties inherent in trying to determine definitive antenna 
positions from attitude data. However, some observations on these 
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results can be made. The procedure can quickly determine the 
system generalized coordinate q^. The coupling of the coordinate 
qg Is quite weak and the study of longer data spans would be re- 
quired to resolve the question of determining qg . The generalized 
coordinates q^ and q^ can probably not be resolved directly by 
use of attitude data alone. Finally, in using the result of these 
runs the best estimate of the true initial conditions would be 
to add the displacement pattern associated with q^ to a symmetric 
equilibrium pattern for q^ . This would yield the following set 
of estimates for the initial in plane displacements: 

A 1 = 122.0 Ag = -158.0 A 3 = -158.0 A^ = 122.0 

4 .3 Cone lusions 

The investigation into the feasibility of determining defini- 
tive antenna positions for the RAE configuration from attitude 
data alone leads to the following conclusions: 

1) It is not, in general, possible to determine 
definitive antenna positions from attitude data 
alone . 

2) The deformations associated with the generalized 
coordinate q^ (Figure 8) are readily determined 
from attitude data. 

3) It may be possible to determine motions of the 
generalized coordinate qg (Figure 8) if they are 
large enough and good quality data exists for 
long periods of time (several orbits). 
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4 ) Small motions associated with the generalized 
coordinates q^ and (Figure 8) cannot be 
determined from attitude data. 
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